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Numerous treatments of the diffuse reflecting properties of scattering media have been 
described. Many theories give an adequate account of the reflectance for a specific set of 
conditions for which the model was constructed and the solution tested experimentally. 
Only those models which are considered to be fairly general are considered here. 

It is convenient to divide the theories into those based upon continuum models and 
those based upon statistical models. The continuum models typically describe the scattering 
and absorbing properties of a given medium in terms of two phenomenological constants. 
These models may all be regarded as varying levels of approximate solution to the general 
equation of radiative transfer. This provides a convenient basis for comparison of the various 
theories. 

The statistical models are based upon a summation of transmittances and reflectances 
from individual layers or particles. Thus, some assumptions must be made about the nature 
of the fundamental units, and the validity of the ultimate result will depend upon how closely 
these assumptions correspond with reality. Only the statistical models lead to expressions 
from which absolute absorptivities and scattering coefficients can be calculated and related 
to the actual particle characteristics. 

The relationship between the various models will be discussed and the features which 
typify the absorptivity and scattering coefficient according to each will be compared and re- 
lated to the available experimental data. This leads to a consideration of the characteristics 
of appropriate model systems and standards. 

Key words: Absolute absorptivities; continuum models; diffuse reflectance; radiative trans- 
fer; reflectance spectra; scattering coefficients; statistical models. 



List of Principal Symbols Used 

(Note : Where a given letter is used in both capital 
and lower case form (e.g., r,R and t,T), the capital 
letter refers to the macroscopic observable and the 
lower case letter to the corresponding variable for an 
individual particle or layer of the material. A bar 
over a given symbol means the average value for 
that variable) . 

a total absorption of a single particle 

(layer); also — ~— (equation (18)) 

Gaussian weighting function (equation 
(33)) 

Rozenberg constants (equation (23)ff.) 

(a 2 -!) 1 ' 2 (equation (21)) 
TJTi polarization of radiant beam (equa- 
tion (24)) 

fraction of radiation reflected (equation 
(75) )_ 

fraction of radiation transmitted (equa- 
tion (76)) 
Rozenberg constant (equation (29)) 



aj 



a>u 



C< 



ft 



H(n) /^-integral of Chandrasekhar (equation 

(35)) 
I radiant intensity 

V component of source function for self- 

radiation (equation (5)) 
j scattering function (equation (2)) 

J source function (equation (4)) 

k absorptivity 

K Kubelka-Munk absorption constant 

(equations (16) and (17)) 
K lf K 2 Attenuation constants in GureviS layer 

model (equations (10) and (11)) 
I particle diameter or layer thickness 

l f "free" path length of Antonov-Romanov- 

sky (equation (58)) 

L Gurevic constant =^K 2 2 —K* 

m mean number of reflections 

M transmittance of a particle for a single 

pass 
M e radiant flux density 

p phase function (equation (2)) 

P n (jtt) Legendre polynomial 

q Rozenberg constant (equation (29)) 

Q Rozenberg multiple reflection constant 

(equation (26)) 
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a single particle 
a single particle 



r total reflectance of 

(layer) 
r simple reflectance of 

(layer) surface 
r* mean external reflection coefficient for 

side scatter 
r e reflectance of a particle (layer) for 

externally incident radiation 
r t reflectance of a particle (layer) for 

internally incident radiation 
R reflectance of a (macroscopic) layer 

R m reflectance of an infinitely thick layer 

s effective path length (equation (58)) 

5 Kubelka-Munk scattering constant 
(equations (16) and (17)) 

S(fx,fjL ) Kozenberg weak scatter constant (equa- 
tion (29)) 

t total transmittance of a single particle 

(layer) 

T transmittance of a (macroscopic) layer 

u radiation emerging from a particle per 

unit solid angle 

u shading factor (equation (70)) 

u d fraction of radiation emerging from a 

particle in a downward direction 

u s fraction of radiation emerging from a 

particle in a sideways direction 

u u fraction of radiation emerging from a 

particle in an upward direction 

x length 

y Johnson multiple reflection factor (equa- 

tion (51)) 

a absorption coefficient 

fi a/a (equation (23)) 

7 Fassler-Stodolski constant (equation (82)) 

8 Fassler-Stodolski constant (equation (82)) 

6 polar angle 

k attenuation eoefScient=a+o- (equation 

(i)) 

X wavelength 

M cos S 

\ij one of the roots of the Legendre poly- 

nominal P n (ju) 

£ effective hole cross section (equation 

(100)) 

p density 

a scattering coefficient (equation (2)) 

r optical thickness (equation (3)) 

azimuthal angle 

<£ e radiant flux 

co albedo of single scatter (equation (30)) 



I. Introduction 

It is now recognized that diffuse reflectance spec- 
troscopy is a very useful companion technique to 
transmission spectroscopy. Not only can it provide 
absorption data in some cases where transmission 
measurements fail, but for many industrial and re- 
search applications, it may in fact be the preferred 
technique. 



Our discussion will be concerned with the behavior 
of radiation within a scattering medium. For sim- 
plicity we assume the scattering centers, which may 
also absorb radiation, are imbedded in a medium 
which neither scatters nor absorbs. The medium 
may usually be taken as air, although there are 
many other cases of interest in which the refractive 
index of the medium is much greater than unity. 
We will not deal specifically with these cases in the 
present review, nor with those processes which alter 
the frequency of the radiation, such as luminescence 
and fluorescence. 

Scattering takes place under a wide variety of 
conditions. One may be concerned with the glowing 
photosphere of the sun, which is surrounded by a 
cloud of electrons that reradiate the direct sunlight 
incident on them. Sunlight is also reradiated by cos- 
mic dust, which accounts for the outer part of the 
corona. 

In a more down-to-earth situation one may be 
dealing in the laboratory with the spectroscopy of a 
powder, paper, opal glass, photographic emulsion, 
etc. As so often happens, concepts developed in one 
area of science are slow to find their way into another. 
We will attempt to show that there is a close rela- 
tionship between the astrophysical solutions, which 
are based largely on radiative transfer theory, and the 
various models which are more familiar to the 
spectroscopist. 

Since we will not be following a historical develop- 
ment, it may be useful to point out some relation- 
ships between the early studies. The first attempt to 
account for transmission and reflection of a layered 
material was carried out by Stokes in about 1860 
[l], 1 and led to some very useful relationships which 
have also been derived by other workers (vide infra). 
Lord Kayleigh [2] and Mie [3] developed the theory 
of single scatter to a high degree, but Schuster [4] 
was the first to consider multiple scatter. He was 
concerned with the cloudy atmospheres of stars, and 
developed a plane-parallel plate model in which the 
radiation field was divided into forward and back- 
ward components. This same model was used much 
later by Kubelka and Munk [5], whose names are 
usually attached go it by spectroscopists. Schwarzs- 
child [6] showed that the radiation field should be 
characterized by a complete angular distribution, 
and if one integrates over the forward and backward 
hemispheres, the Schuster model is obtained as a 
first approximation. A further generalization of the 
Schwarzschild formulation leads to an integro- 
differential equation known as the equation of radia- 
tive transfer, which is very general in concept, but 
can be solved exactly in only a few cases. 

The radiative transfer theory and various models 
stemming from it are referred to as Continuum 
Models. They have in common the characterization 
of the scattering and absorbing properties of the 
medium through phenomenological constants, usually 
two in number. These theories will be con- 
sidered in section II. A completely satisfactory theory 
must of course relate the measurable quantities to 



1 Figures in brackets indicate the literature references at the end of this paper. 
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fundamental properties of the scattering material, 
such as particle size, refractive index, and absorpti- 
vity. This is the goal of the statistical theories. 
Many have been proposed; those which show the 
most promise are discussed in section III. In section 
IV we discuss some of the strengths and weaknesses 
of the various models and some relationships be- 
tween them. Section V considers the meaning of the 
absorption and scattering coefficients often used to 
interpret reflectance data, as well as the character- 
istics required of appropriate model systems. 

We will limit our discussion to ideal, homogeneous 
dispersions. There have been interesting develop- 
ments in the theory of nonhomogeneous media, mix- 
tures, luminescing materials, diffusing media in 
which photochemical reactions are taking place, and 
reflectance for materials dispersed in a highly re- 
fracting matrix, but these topics are considered to 
be outside the scope of the present review. Many 
of these topics are covered in the books which have 
been written on the subject [7,8,9]. Further informa- 
tion has been presented at various symposia which 
have been held, and the proceedings of at least two 
of them are available [10,11]. These sources, together 
with various reviews devoted to applications in 
particular fields, have made workers aware of the 
power of the technique. Thus we will confine our- 
selves to certain theoretical aspects. 

The present paper may be regarded as an update 
to a former review of the subject by the author [12]. 
Most of the theoretical work of the intervening years 
has been in the development of more refined statis- 
tical models, and in showing the relationship be- 
tween the various theories. At present there is some 
reason to feel optimistic that out of a morass of 
apparently divergent and unrelated theory, some 
order is beginning to emerge. 



II. Continuum Theory 

As is often the case with applications of physical 
theory, the real media with which diffuse reflectance 
spectroscopy is concerned are intermediate between 
two extreme cases, each of which is well understood. 
The one limiting case involves the propagation of 
radiation through quasi-homogeneous matter, where 
the radiation field is characterized by smoothly and 
slowly varying functions of the coordinates. The 
other extreme involves the emission, absorption, and 
scattering of radiation by single particles in a 
homogeneous medium. The local transformation of 
the radiation field at the boundaries of quasi- 
homogeneous media is a special case which has 
likewise been characterized for a long time, and leads 
to such well-known phenomena as reflection, refrac- 
tion, etc. 

In dense scattering media, it is important to 
consider the manner in which the radiation fields 
from the various scattering centers interact with one 
another. Rozenberg [13, 14, 15] has pointed out that 
the interaction can be treated as a sum of two parts, 



one of which is coherent and the other incoherent. 
The coherent part is largely due to nearest neighbor 
interactions and gives rise to dispersion effects; i.e., 
is involved with changes of refractive index. The 
incoherent part is that with which we are primarily 
concerned in the present discussion. It involves 
multiple scattering as a sum of interactions through- 
out the entire medium. In most treatments the 
multiple scattering, because of its complexity, is 
separated rather arbitrarily from dispersive effects. 
Such a separation reduces the problem to one involv- 
ing geometrical optics, and allows one to write a 
radiative transfer equation in which the absorbing 
and scattering properties of the medium are treated 
as phenomenological constants (spoken of as "ex- 
ternal parameters" by Stepanov[16]). 

The concept that each volume element is irradiated 
by scattering from every other volume element of 
the medium (The Principle of Self-Illumination), is a 
basic concept of radiative transfer theory which 
clearly pervades the writings of such early workers 
as Schuster [4, 17], King [18], and Schwarzschild [6]. 
The equation of radiative transfer can in fact be 
regarded as a simple statement of the law of con- 
servation of energy. The change in intensity of a 
beam along its direction of propagation, dl, is equal 
to the radiation which is lost through absorption and 
scattering —Kpldx, plus that which is scattered in 
this direction from all other directions, jpdx: 



dl= —KpIdx-\-jpdx 



a) 



Here p is the density, k is the attenuation coefficient, 
j is the scattering function, and dx is the element of 
path length. The scattering function can be written 



as 



M *)=;£ £ f**P( e > # *'> 4>')W, *') sin O'dB'W 

(2) 

where a is the scattering coefficient and p(d,<t>; 
6' ',$') is the phase function which defines the 
probability that radiation which travels initially in 
the direction (0',</>') is scattered into the direction 
(d,^). Equation (1) is usually solved in terms of 
the optical thickness 



rx 2 
t= Kpdx 

Jxi 



(3) 



Defining as the angle with respect to the inward 
surface normal gives 

where /x = cos and J=j/<r is the source function. 

The equation of radiative transfer can be general- 
ized to include dispersion effects as well as scattering 
[14, 15]. In terms of the components of the Stokes 
vector, I, 
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dx 



i-i v. 



0) 



+|" I ^Pt 6 ' * ; B '> *') 7 ^'> *') sin o'M'W \ 

+/«'(»,*). (5) 

The term Z"/ (0,0) accounts for possible self radiation 
of the volume element, which is of importance in 
the infrared spectral region or in luminescing media. 
As mentioned above, this term will be suppressed in 
the present discussion. 

A. The Layer Model 

A model in which the scattering medium is ap- 
proximated by plane-parallel layers has been used 
by numerous workers as a basis for reflectance 
theory [7, 8]. We consider a pair of adjacent layers, 
with RijRj and T u Tj as the reflectances and trans- 
mittances for radiation incident in one direction, 
and RijRj and T I} Tj as the reflectances and trans- 
mittances for radiation incident in the other 
direction. 

If the incident beam is of unit intensity, then the 
portion R t is reflected and the portion T t transmitted 
by the first layer. At the second layer the portion 
T t Rj is reflected and T t Tj is transmitted. The beam 
TJZj strikes the underside of layer i where TJtjTj 
is transmitted, while TiRjRj is reflected. Continuing 
this process indefinitely (see fig. 1) we find that the 
transmittance and reflectance of the combined 
layers are given by 

T i+j =T i T j (l+R I R j +R I 2 R j 2 + . . . ) 

R i+j =R i +T i T I R j (l+R I R j +R I 2 R j 2 + . . . ) 

Summing the geometric series gives 



Ti+j— y_ 



tii+j — Ri 



-RjRj 

■ TjTjRj 
1 — RjRj 



(6) 



(7) 
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Figure 1. Reflectance and transmittance of a pair of inhomo- 
geneous layers (Kubelka [19]). 



We generally deal with the case where Ti=Ti 
and R i =R I ; i.e., where the reflectance and trans- 
mittance of a layer do not depend on the direction 
of incidence. This is quite a general condition 
which applies even to the case of nonhomogeneous 
sheets [19]. With these substitutions eqs (6) and (7) 
reduce to 

T T 

(8) 



Ti+j—-^_ 



-RiRj 



Ri+j—Rc 



T t *R, 
1 — RiRj 



(9) 



We now make the assumption that the layers are 
homogeneous and thick compared with individual 
particle dimensions so that no intrinsic inhomogene- 
ities are evident. We can then write that the reflec- 
tance of a given layer is proportional to its thick- 
ness [20], 

R(dx)=K l dx. (10) 



Similarly we write 



T(dx)=l-K 2 dx. 



(11) 



Here K 1 and K 2 are phenomenological constants 
which characterize the medium. We assume that 
K 2 >K U where the inequality applies to absorbing 
media and the equality to nonabsorbing media. 
We substitute R(x) and T(x) for R t and T t in eqs (8) 
and (9), with R(dx) and T(dx) from eqs (10) and (11) 
being substituted for Rj and Tj. Writing R i+j and 
T i+j as R(x + dx) and T(x+dx) allows eqs (8) and (9) 
to be expanded in series to give 



and 



dR=K 1 T 2 dx 
dT=-T(K 2 -K 1 R)dx. 



Integrating this set of coupled differential equa- 
tions subject to the boundary conditions 2?(0)=0, 
T(0) =1 gives 



R=R n 



l—e~ 2Lx 
\-Rte~ 2Lx 



T=(l-Rl) 



(12) 



(13) 



R w is the reflectance of an infinitely thick layer, and 
it is related to Ki and K 2 through 



R, 



K 2 -L Ki-y/Kf-Ki 2 



K x 



K, 



while L is given by 



Z=V2f 2 2 -^i 2 



(14) 



(15) 



Equations (12) and (13) were derived long ago by 
Stokes [1] for plane-parallel plates, with similar sub- 
sequent derivations by Schuster [4, 17] and Gurevic 
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[20] for light-scattering layers. These are now known 
to be a special case of the more general Kubelka- 
Munk theory [21], which we will now consider. 



B. The Kubelka-Munk Theory 

The Kubelka-Munk theory [5] is based upon a 
model in which the radiation field is approximated by 
two fluxes, the one, i+, traveling from the illuminated 
sample surface, and the other, /_, traveling toward 
the illuminated surface (see fig. 2). As radiation 
travels from the surface, its intensity is decreased by 
scattering and absorption processes, both assumed 
to be proportional to the thickness of the medium 
traversed. This is partially offset by scattering from 
the other beam, so we have 



dh 



-(S+K)I + dx+SI_dx. 



(16) 



The component travelling toward the illuminated 
surface is similarly described: 



dI_=(S+K)I_dx-SI+dx. 



(17) 



R- 



_(R t -R a )/R a -R a (R-l/RJ exp [SX(l/R„-R a )] 



R-R a -(R-l/RJ exp [SX(l/R a -RJ] 

(20) 

R g is the reflectance of the background (see fig. 2) 
and R a is once again the reflectance of a layer which 
is so thick that further increase in thickness does 
not alter the reflectance. If eq (19) is integrated over 
the limits x=0 to x=™ , a simple formula results [5] 



2? =limi2= 



.(a 2 -l) 1/2 =a-b 



(21) 



Using eq (18) this can be rearranged to give the 
well-known Kubelka-Munk function F(R m ) 



*Bj.fl=5JLA 



2R„ 



S 



(22) 



We note in passing that this result follows directly 
from eq (12) in the limit x— >°°[23], which once again 
shows the close relationship of the Gurevic and 
Kubelka-Munk models. 



The constants which we have introduced here are 
once again phenomenological constants which 
describe scattering (S) and absorption (K) within 
the medium. If we make the following definition, 



a = - 



S+K 
S 



(18) 



we can write 



-dl 4 



Sdx 

dl- 

Sdx' 



—al++l. 



which can be combined into a single differential 
equation, 

M-&-20B+1 (19) 

where R=IJI + . Equation (19) can be easily inte- 
grated over the entire thickness x of the scattering 
medium to obtain [5] 



v '■'.'• 1 1- 
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dx 



Figure 2. Model for the Kubelka-Munk analysis of reflectance 
and transmittance of a scattering medium [5]. 



C. Rozenberg Solutions 

For a homogeneous semi-infinite medium, a good 
approximate treatment of the reflectance for strong 
absorbers has been given by Rozenberg [13-15]. The 
solution involves summing various successive con- 
tr butions to the reflecting power by scattering of 
different degrees of multiplicity, and is based upon 
concepts developed by Kuznetsov for problems of 
visibility [23]. To the nth degree of approximation, 
the reflected intensity in an isotropic medium is given 
by 



«,*)=/,•§!!#• 



(23) 



In this case 6 and <t> define the direction of obser- 
vation, 1 1 is the incident beam intensity, 0=— is the 

ratio of the absorption and scattering coefficients, 
and the index i runs over the four components of 
the radiation field (see eq (5)). The a u coefficients in 
eq (23) are given by 

a«=S*w(M)0?(*o,*>) 

where the C5(0 o ,*o) = *?/*! denne the polarization 
of the incident beam, and the a ijt (0, 0) are coeffi- 
cients which depend only on the angles of incidence 
(0 o ,0o) and observation (0,0), and on the form of the 
scattering indicatrix. 

The reflectance of the medium is given by 



C t R= 






(24) 
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where Ct=Ii/Ii is the polarization of the reflected 
radiation. It will be observed that as increases, 
scattering of higher multiplicities becomes less im- 
portant. Ambartsumian [24] has shown that the 
mean multiplicity of scatter in the case of reflection 
from a semi-infinite turbid medium is Vl + 1/0. Thus 
when > 1/3, a fairly accurate solution is obtained by 
inclusion of terms up to second or third degree. It is 
further assumed that the scatter is independent of 0, 
which should be a good approximation for mixtures 
of polydispersive media with different a and a; i.e., 
as with the addition of dye to a suspension. The re- 
sulting equation is 



flo ji+0) 2 
R i+i- 



(25) 



where R is the reflectance of the medium itself 
(when 0=0 O ), and 



y 1+ &an(l+A))' 



(26) 



Thus Q is a quantity which defines the relative 
contribution of higher multiplicities of scattering 
when j3=0 o . Both it and R are constants which are 
independent of the nature and concentration of the 
colorant. 

Equation (25) may be regarded a^ a generalization 
of Lambert's law to the case of colored media. It 
has been derived in a somewhat different form by 
Chekalinskaia [25] from scattering theory. In terms 
of the reflection (r), forward scatter (t), and absorp- 
tion (a) constants of a single scattering layer( a-\-r 
+£=1) used by Chekalinskaia, the Rozenberg con- 
stants can be written 



0: 



a 



<2=14 



t 



r+t 



(27) 



D. Exact Solutions 

In problems of spectroscopy we often assume 
isotropic scatter. We know that in no case is single 
scatter actually isotropic [3], although the random 
distribution of anisotropic particles and scatter 
apparently tends toward an isotropic result [15]. 
Problems in highly anisotropic scattering media 
have been considered by some workers [29-34]. A 
detailed discussion of these solutions, often obtained 
by numerical computer methods, will not be discussed 
here. 

For the simple case of isotropic scatter, the 
phase function (see eq (2)) can be written 



p(M;«'tf')=«o= 



<j-\~ol 



(30) 



Here co is known as the albedo of single scatter. 
It represents the fraction of the radiation lost by 
scattering in a medium where both absorption (a) 
and scattering (a) take place. With this assumption, 
the equation of radiative transfer for a plane- 
parallel semi-infinite medium becomes independent 
of the azimuthal angle $ and we have 



dr 



r+i 



=^(t,m)— § w °| I(T )f j.')dn' . 



(31) 



The integral occurring in eq (31) may be approxi- 
mated by a Gaussian quadrature, in which case a 
set of coupled linear differential equations is obtained, 



Vi -T^=l i—7i ^o S Qjlj- 



di 



S 

j = — n 



(32) 



The constants a ; are Gaussian weighting functions 
given by 

1 C +1 Pn(l*)j /ooN 

a j= V it \ d V \oo) 

r n KM J -l At" fij 



Il'ina and Rozenberg [26] have demonstrated the 
validity of eq (25) in several instances. Obviously 
in a highly absorbing medium where 0>>1, eq (25) 
further simplifies to 



t-<* 



(28) 



In the other limit; i.e., where 0<<1, Rozenberg 
[15, 27] has shown that the reflectance can be written 
as an exponential of the form, 

tfco U Mo)=— Mm, mo) exp [-4Vg0S( M , Mo)] (29) 



where h(fjL, /z ), S (/z, Mo)> and q are quantities which 
depend on the form of the scattering indicatrix. Ro- 
manova [28] has determined these quantities by exact 
solution of the radiative transfer equation. 



and iJij is one of the zeros of the Legendre polynomial, 

P»(m). 

Passing to the limit ?i— >oo gives an exact solution 
which Chandrasekhar [35] has shown to be of the 
form, 



J r (M)=j°*.U)^ffG.)ffU). 



(34) 



Here I r {y) is the reflected intensity in the direction 
fi from a collimated incident beam in the direction 
ix , whose flux per unit area normal to the beam is 
ir$ e (iM,)- The H— integrals are denned by 



ff(M)=l- 






r 



n+n 



(35) 



and tables of them have been given by Chandrase- 
khar [35]. 
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Giovanelli [36] has given explicit expressions for 
several cases of interest. The total reflectance for 
light incident in the direction /x is 

i?U) = l-i?(Mo)(l-co ) 1/2 (36) 

while that for diffused incident radiation is 

R D =l-2(l-u ) 1/2 (\Hb)dv. (37) 

Jo 

Tables of the first moment of the H — integral, 
which occurs in eq (37), have also been given by 
Chandrasekhar [35]. 

Exact solutions to the equation of radiative 
transfer can be derived for other phase functions 
as well. In general, the phase function may be 
expanded as a series of Legendre polynomials 



p(cos Q)=^DcoiPi(cos 0) 

1=0 

where axial symmetry is assumed. Terms higher 
than first degree contribute very little [36, 37], 
and thus the approximate phase function. 

p(cos e) = co (l+x cos 0)(O^z^l) 

is sometimes used. Exact solutions are available 
for scatter according to this phase function also 
[36]. 

Equations (36) and (37) for isotropic scatter can 
be readily applied using tables given by Giovanelli 
[36], and they of course are of considerable theoretical 
interest since they represent exact solutions to 
which the various approximate theories can be 
compared. We do not expect all media to scatter 
isotropically, but we might expect the range of 
applicability of the equations to be extended if an 
appropriate average scattering coefficient were used. 
If we consider a diffuser in which each scattering 
center scatters light symmetrically about the direc- 
tion of incidence, we may write [38] 



where 



<Teff=(l — mV 



r+i 



(38) 



(39) 



In this approach it is assumed that the same isotropic 
solutions (eqs (36) and (37)) may be used for 
arbitrary angular distributions of scatter, so long as 
the scattering is averaged according to eqs (38) and 
(39). Blevin and Brown [38] have shown that the 
reflectance curves are essentially the same for 
isotropic scatter or for scatter according to the 
phase functions l+Pi(/z), 1+^Mm), an d l+PaGz). 
This suggests that the reflectance is not a sensitive 
function of the scattering indicatrix, and the isotropic 



solutions are in fact a good approximation for res 
scattering media. 



III. Statistical Theory 

Continuum models, as we have seen, are somewhat 
limited. They involve the use of phenomenological 
constants with no obvious relationship in general to 
the fundamental constants with which we are familiar 
(molar absorptivity, refractive index, particle size 
and shape, etc.). Statistical theories, on the other 
hand, involve the construction of an appropriate 
model and the success of the theory depends on just 
how closely the model approximates real sample 
conditions. 

It appears certain that one of the most severe 
limitations of continuum models is the assumption 
that they remain valid even when infinitesimal thick- 
nesses are considered. This is in fact contrary to the 
assumption of homogeneous layers previously in- 
voked (see section II-A) , and it is this contradiction 
which is largely responsible for limiting the range of 
applicability of continuum models, as our subse- 
quent discussion will show. 

Let us return to eqs (8) and (9) and assume that 
we are now dealing with thin layers whose thickness 
is that of the individual particles. If we take layer i 
to be the first layer and layer j to be the combination 
of all the other layers of an n-layer sample, we have 



k 



ht 



2,3, ... n~ 



2, 3, 4, ... n 



1 rir 2f3>4f 



and 



7*1, 2, 3, . . . n — r l~ 



ti ^2, a 



1 rir 2| 3, 4, 



(40) 



(41) 



Passing to the limit n-+ oo } we write 

^1, 2.3. . . . n~ ^2,3.4, . . . n~ * co— 
r l, 2, 3, . . . n ==r 2, 3, 4, . . . n = R<o 

Equation (41) then becomes 



R m =r- 



l-rR a 



(42) 



where we have assumed that all layers are the same 
so the subscripts on r and t can be dropped. Equation 
(42) can be solved for R m to give an expression for 
the reflectance of an infinitely thick sample in terms 
of the reflectance and transmittance of a single layer. 
The result is 



R a 



l+r 2 -t 2 
2r 



-V(^7- <« 



This equation is fundamental to essentially all sta- 
tistical theories, the only difference being in the 
method used to calculate r and t. 
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We have seen that the Kubelka-Munk 
leads to a solution of the form, 

(i-R m y__K 



2B a 



"8 



When this is solved for R m , we get 



R a 



_K+S-^(K+S) 2 -S 2 
S 



theory plot of F(R m ) versus K deviates from linearity for 
high values of K [7—9], and it appears that eq (48) 
can be used to explain the deviations in part. It 
should be recognized that the deviations at high 
values of K are probably a result of anomalous dis- 
persion effects also, but eq (48) does represent an 
improvement in the range of validity and shows the 
need to consider the particulate nature of scattering 
media in developing a more precise theory by which 
(44) absolute absorptivites can be determined. 



which is not of the same form as eq (43). 

We assume with Simmons [39] that the plane- 
parallel layers of the Kubelka-Munk model cannot 
be made infinitesimally small, but are restricted to 
layers of finite thickness I, where I may be inter- 
preted as the mean particle diameter of the sample. 
Then the fundamental differential equations of the 
Kubelka-Munk theory (eqs (16) and (17)) are re- 
placed by the finite difference equations: 



di + ^(i + ) t+l -d + ) i 



dx 



I 



= -(K+S)(I+) i +S(I.) t+1 (45) 



<*/__(/_)„.,-(/-)< 



dx 



/ 



= (K+S) (I-) i+1 -S(I + ) ( (46) 



where the subscripts i and (i+1) refer to the ith 
and (i+l)st sample layers, respectively. Now for an 
infinitely thick sample, 

*• (/+)< (i + ) t+1 

and Eqs (45) and (46) can be solved to give 



Rm = 

2(S+K-KlS-KH/2) 



-^(S+K-KlS-KH/2) 2 -AS 2 



2S 



(47) 



Equations (43) and (47) are identical if we make the 
following identifications : 

S=r/l 

K=(l-r-t)/l=a/l 

where a is the fraction of the incident radiation which 
is absorbed by the layer. 

The difference between the traditional and modi- 
fied Kubelka-Munk solutions may be seen by writing 
eq (47) in the form, 



K K 2 1 

F(R„)=^-K-^=a 



\r J 2r 



(48) 



It will be recognized that the difference lies in the 
addition of the last two terms. It is well known that a 



A. The Bodo Model 

Bodo [40] used a procedure similar to that used to 
derive eq (43) for the derivation of r and t. We will 
denote the simple reflectance of the layer surface 
by r , the absorptivity (defined through I=I exp 
(—kx)) by k, and the layer thickness (equivalent to 
the mean particle diameter) by I. Then according to 
figure 3, the reflectance and transmittance of a single 
layer are given by 

r=r + (1 -r ) 2 d- 2kl +(l-r ) 2 r z e- ikl 

+ (l-r ) 2 r Q 5 e-« kl + . .. 

t=(l—r )%- kl +(l—r )\ 2 e- m +(l-r )We~ m 

+ (l-r ) 2 r 6 6- 7 **+ .... 



Summing these series gives 

_r [l + (l-2r ) exp (-2kl)} 



t= 



1— r 2 exp (-2kl) 

(l-r ) 2 exp (-kl) 
1— r 2 exp (—2kl) 



(49) 



(50) 



Equations (49) and (50) together with eq (43) 
constitute the Bodo formulation, which is in fact 
equivalent to that of Stokes [1] and Girin and 
Stepanov [41]. Bodo [40] obtained good results with 
these formulae for powdered glass samples using the 
arbitrary assumption that r =0.10. Karvaly [42] has 
shown that this was at least in part due to a particu- 
larly favorable position for the absorption band 



/, ,\2 -2lcl 

(1 - r ) 'o e 



<1-r ) 2 r.V 4kl 



O-OVe" 6 " 




. .2 2 -3kl 

(1- r ) r e 

,, .2 4 -5kl 
(Wo) r o e 



. .2 6 -7kl 

(1 - r ) r e 



Figure 3. Reflectance and transmittance of a single layer of 
thickness I according to Bodo [40]. 
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chosen for study, but in the general case, R m is a very 
sensitive function of r . 

Bauer [43] showed that in some cases the layers 
should be considered to have rough surfaces where 
total internal reflection can take place, and he has 
derived expressions analogous to eqs (49) and (50) 
for this case. 

B. The Johnson Model 

Johnson [44] has carried out the summation some- 
what differently than Bodo, but with quite similar 
results (see fig 4). It is assumed that there are p 
layers and that the mean number of attenuating 
reflections which the rays undergo in the 2p traver- 
sal is yp, so that the reflectance is given by 

R p =r +2r J2(l-ro) vp exp (-2kpl). (51) 
V 

Thus y can be regarded as an adjustable parameter 
which gives a semi-empirical account of multiple 
reflections as well as scattering losses. The sum for 
an infinite number of layers is 



B.=r +2r 



exp [y In (1— r )— 2kl] 



1— exp [?/in (l—r )—2kl] 
which is equivalent to eq (43) with 

r=r [l+(l-r )*exp(-2K)] 
Ml-r ) v/2 exp(-JfcZ). 



(52) 

(53) 

(54) 



The denominator of equation (52) can be expanded 
to srive 



fi« 



=r 



'2(l-r ) y exp (-2H) 



p(i-r : 

L 2kl- 



2/ In (1— r ) 



a 



(55) 



This result could also be obtained directly by inte- 
gration of eq (51), which suggests that it may in fact 
be a more realistic representation of a real sample 
whose particles actually have a range of diameters. 



\2 -2kl 

l*o(l -'d) e 



)2r (l-r ) 3 *e- 6kl 




It will be observed from figure 4 that y should be 
set equal to 4 for the case of no multiple reflections. 
In such a case only one half of the incident light is 
reflected, however. Johnson has suggested that y can 
be estimated from eq (55) by setting R co = l for 
k=0. This gives 



(l-ro)'- 1 ! , 1 1 , 1 3 , 

y - 2 + 4 ro, '6 ro+ 8 r ° + - 



(56) 



from which y=2 is seen to be a satisfactory approxi- 
mation for refractive indices smaller than 1.5. A 
smaller value of y is required for larger refractive 
indices. Companion and Winslow [43] have used a 
model similar to Johnson's, but which includes all 
multiple reflections. The summation was carried out 
by computer and no explicit expression for the re- 
flectance was given by these workers. 

Johnson [44] also suggested that r be equated to 
1.5 times the normal Fresnel reflectance. This is 
meant to account for the random distribution of 
particle surfaces and corresponds with an average 
incidence angle of approximately 30°. It was shown 
[44] that eq (55) yields absorption coefficients for 
KC1:T1, KBr:Tl, and didymium glass which agree 
satisfactorily with those obtained by transmission 
measurements of the same materials. 



C. The Antonov-Romanovsky Model 

Antonov-Romanovsky [22] has developed expres- 
sions which can be used to calculate the true absorp- 
tion coefficient from reflectance measurements by 
connecting the Kubelka-Munk and Bodo theories. 
Antonov-Romanovsky treats two limiting cases of 
regularly-shaped sample particles, spheres and paral- 
lelepipeds (see fig. 5). 

For spherical particles the radiation impinges on 
the surface from within at the same angle that it 
entered the particle, since the angles which the cord 
of a circle makes with respect to the surface normal 
must be the same. Therefore, total internal reflection 
is impossible, and I will approximate the particle 
diameter, I, actually being somewhat smaller. 

In the case of the parallelepiped, some total 
internal reflection is possible, but most radiation 
probably exits through an opposite face without 
further reflection. In this case also it is obvious that I 





Figure 4. Model for Johnson analysis of reflectance from a 
scattering medium [44]- 



Figure 5. The Antonov-Romanovsky model for regular spheres 
and parallelepipeds [22]. 
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is approximately of the same magnitude as I, but 
somewhat larger. 

Thus for regularly shaped particles it is assumed 

that Z~ I and 

, 1, (l-r ) 2 -2r (l-2r )^(iQ 

^~2Z ln (l-r y-2r F(BJ {b7) 

This is an approximate form of the Bodo model which 
is valid for media in which r 2 <C<O- 

For irregularly shaped particles it is assumed that 
the emerging radiation meets the surface of the 
particle with an equal probability for all angles. This 
will be the case if the mean number of reflections in 
the layer m , is large : 

m>l (7&>1.5). 

This condition requires that 

where l f is defined as the "free" path length. The 
effective path length is then 

s=rhl f (58) 

and it is assumed by Antonov-Romanovsky that 

l f ~l/2. (59) 

The assumption that the emerging radiation is 
independent of angle allows us to divide it equally 
between the two sides of the layer whose thickness is 
I, the mean particle diameter. We write 



r— r =t 
and from the law of conservation of energy, 

(l-r )(l-exp(-sk)) = l-r-t. 
Equations (60) and (61) can be solved to give 
_2r +(l— r ) exp ( — mkl/2) 



r =- 



(1— r ) exp {—mkl/2) 



Using these in eq (43) gives 



2_ (i-r )'-(l-r,)F(fl,) 
ml m (l-r y-2r F(R„) 



(60) 
(61) 

(62) 
(63) 

(64) 



D. The Melamed Model 

^ In contrast with the statistical theories previously 
discussed, Melamed carried out a summation over 
the reflectance and transmittance of individual 
particles rather than layers [46]. Some features of 
this model are better understood using an alterna- 



tive derivation of Karvaly [47], which we follow 
here (see fig. 6). A single particle of the surface layer 
is shown in the figure; it is shifted laterally with each 
reflection to illustrate the path of the radiation. 

t Of the diffuse radiation which is incident on a 
given surface particle, the fraction 2ur e is reflected, 
where the subscript e is used to indicate reflection 
of externally incident radiation. Here u is defined as 
the radiation emerging per unit solid angle from a 
particle. 47ru is the solid angle which would be 
observed from the second layer if a particle were 
removed (see fig. 7). 

Of the radiant flux (l-2w e ) which enters the 
particle, the part (l-2ur e )ut contributes to the re- 
flectance from_the first internal reflection, and the 
part^L = (l-2ur e ) (l-u)t strikes the underlying layers, 
which are considered to form an infinitely thick 
powder mass of reflectance R a . The transmission 
of an individual particle is represented by the 
symbol t. The part A R m is reflected back into the 
layer,_ where A R m (l-r e ) enters the particle and 
A R m r e is reflected back into the underlying layers. 
Of that part which entered the particle, AoR^il- 
r_ e )ut contributes to the reflectance and AJR m 0-- 
r e )(l-u)t is reflected downward again. This combines 
with that externally reflected to give A R m r e + 
A R m (l-r e )(l-u)t=A R (a [7 e +(l-u)t]^A R a> Q as the 
flux reflected downward of which A Q RJQ returns 
to the layer, etc. 



\ t / 



/ 



/ 



1 2ur e (l-2ur e )ut A RJl-r e )ut AoR.Qd r e )ut 




Figure 6. The Melamed Model for powder reflectance as 

viewed by Karvaly [47], 




Figure 7. Illustration of the meaning of the radiant intensity 
factor u [47]. 
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Proceeding in this manner we can write the follow- 
ing expressions for the reflectance and transmittance 
of a layer for radiation from outside (r h ti) and 
inside (77, t T ) the sample: 

r 1 =2ur e +(l—2ur e )u4 (65) 

r 7 =r i +(l-r i )(l-u)« (66) 

t 1 =(l-2v? c )(l-u)t (67) 

t T = (l-r e )u.t (68) 

Putting these expressions into (cf. equation (7)) 

for an infinite number of layers gives 

K oa =2ur<,+ 7 - — p - N — yz — -tjz — - s P , (69) 
(1 —li m r e ) — (l—u)(l — r e )E m t 

The application of eq (69) requires some assump- 
tions. Melamed assumes the sample to be composed 
of randomly-shaped particles which can be approxi- 
mated as spheres. Then if we assume the absorbance 
is small (kl<l) and that the radiation leaves the 
sphere isotropically, 



u= 



U 



l — (l — 2u )t 



(70) 



where the shading factor u for close-packed spherical 
particles is 0.284. The transmittance of a single 
spherical particle including multiple internal reflec- 
tions is written in terms of the internal reflectivity r f : 



t= 



(l-r,)M 



1— tM 



where 



M= 



(kiy 



■[l-(fcZ+l)exp(-4Z)] 



(71) 



(72) 



is the transmittance of the particle for a single 
pass^ [48]. The mean external and internal reflec- 
tivities, r e and r u were calculated by Melamed as 
the averages over the Fresnel reflectivity r(0), 

rr/2 

r e =2 r(6) sin 6 cos Odd (73) 

r<=(l— sin 2 6 c )+2 r(0) sin 6 cos Odd (74) 

for an ideal diffusing surface obeying the Lambert 
cosine law [49]. Here 6 C is the critical angle. 



E. The Fctssler-Stodolski Model 

Fassler and Stodolski [50-52] have pointed out 
that the Melamed theory is not only cumbersome to 
apply, but contains an inconsistency as well. It is 
assumed in the Melamed treatment that the radia- 
tion is distributed isotropically. However, we have 
seen in the previous section that when the Melamed 
equation is expressed in terms of the general layer 
model, we require different reflectances and trans- 
mittances for external and internal radiation (see 
eqs (65-68)). 

A model which preserves some of the features of 
the Melamed theory, but which removes the incon- 
sistency, has been constructed by Fassler and Stodol- 
ski as follows: we assume radiation to _be externally 
incident on a particle layer. The part r e is reflected 
at the surface, and the part (1— r e ) penetrates the 
layer. If we use the symbol a to represent the part 
which is absorbed in the layer, including the effect 
of multiple reflections, then the part (l—a)(l — r e ) 
leaves the particle either above or below the layer. 
We use j r and j t to represent the fractions of the 
radiation directed upward and downward, where 

Then the reflectance and transmittance of the layer 
are given by 



r=r,+/ f (l-o) (1-r.) 
t=f t (l-o) (1-r.) 



(75) 
(76) 



For weakly absorbing systems it is assumed that 
fr=ft, corresponding with the Antonov-Romanovsky 
treatment of irregular particles. In the general case, 
/r</«. Equations (75) and (76) with/ r =/,=K can be 
put into eq (42) to get 



**■>-<■-»•>' i+v-a-?.). 



(77) 



F(R) a is once again the Kubelka-Munk function. 
Solving eq (77) for a gives 



F(B m ) 



(1+7 e) 

"(1-r.) [(l-r.)+F{R m )] 



(78) 



If we assume the radiation within the powder 
layer to be isotropic, we can write (l—t) as the 
fraction absorbed, where t represents the trans- 
mittance of a given particle taking multiple reflec- 
tions into account. Letting u s represent the fraction 
of the radiation from a given particle which remains 
within the layer, we have u s t as the part which 
enters a neighboring particle where (l—t)u s t is 
absorbed. Continuing in this way we once again 
generate a geometric series which can be summed to 
give the absorbance of the layer': 



a= 



l-t 
l—u s t 



(79) 
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Following Melamed [46] the radiation is distributed 
among the fractions going up (u u ), going down (u d ), 
and going sideways (u s ) : 

u u +u d +u s =l. 

Since we assume an isotropic distribution of the 
radiation, we can write u d =u u =u , and so 



u s 



= 1— 2Un 



where once again u =0.284 for spherical particles. 
Using this in eq (79) gives 



a= 



(H 



(i-l+2».) 

Introducing eq (71) we can express this result in 
terms of the transmittance for a single pass, 



a=- 



ir> 



(80) 



^-1+2^(1-?,) 



The Duyckaerts equation for M (eq (72)) leads to 
a very complex expression, but with the Felder 
approximation [53], 

M=ex P r-|Hl (81) 

equations (78) and (80) can be combined to give 

(82) 



where 



and 



M~2 [n l-8F(R„) 



_ 2u (l-r i )(l + r e ) 

7 (1-TeY 



8= 



27, 



(l-r e ) 2 



The parameters 7 and 8 may be treated as 
constants which depend only on the refractive index, 
n. WhenH is small, (y-8)F(R a )<z:i and 8F(R m )<^l 1 
and eq (82) simplifies to 



kl=±yF(R„). 



(83) 



Thus there is a simple linear relationship between the 
true bulk absorption coefficient and the Kubelka- 
Munk function for weakly absorbing materials. 

F. The Simmons Model 

Simmons [54] has used a simplified particle model 
to relate diffuse reflectance to fundamental optical 
constants without the use of the cumbersome equa- 



tions which result from the more refined Melamed 
theory [46]. 

The law of conservation of energy requires that 
the total radiant flux impinging on a particle ($ e ) 
must equal that returning from it ($/) except for the 
part which is absorbed: 



If the spherical particles are assumed to scatter 
isotropically then $>' e = wl 2 M e , where M e is the 
radiant flux density, and 



$ e =Tl 2 M e /(l-a). 



(84) 



Since the reflectance of an infinitely thick sample is 
independent of depth, we can write the two alterna- 
tive expressions for the reflectance, 



R= 



QrPMe/2) 



R~- 



(ttP/2) r l2 M e (6) sin Odd 
(ttZ 2 /2) ( V M e (6) sin Odd 

Jjr/2 

(irl 2 M e /2) 



(85) 



(86) 



<£ e is the total flux integrated over the whole sphere, 
so eq (85) and (86) combine to give 

<S> e =(7rZ 2 /2) {*M e {6) sin 6de=irl 2 M e (R+R- l )/2. 

Equating this to eq (84) gives 

R=[l-(2a-a 2 ) l ' 2 ]f(l-a). (87) 

We let M equal the absorption from a single pass 
through the spherical particle. Of the incident 
radiation the fraction (1— r e ) penetrates the particle 
where the fraction (1— r e ) (1 — M) is transmitted, and 
of this, the fraction 7\- is internally reflected. The 
total fraction absorbed following an infinite number 
of such inter-reflections is 

a=(l-7 e )[(l-M)+r i M(l-M) 

+r/M 2 (l-M)+ . . .]. 
_ (l-r e )(l-M) 
1-rM 

Again using the Felder approximation [53] (equation 
(81)) for spherical particles gives 



(l-r e )[l-exp (-2AZ/3)] 
l-r^exp (— 2H/3) 



(88) 



A simple expression which is valid for small values 
of a can be derived by noting that eq (87) represents 
the leading terms of a series expansion of the form, 



R=l-(2a) 1/2 +a- 



=exp [-(2a) 112 ]. (89) 
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If the exponential terms in Equation (88) are likewise 
expanded we find 

(l-r e )2H/3 ^ l-T e 2kl 



a~-^= 



l—r t +2klr i /3 1— r< 3 



It can be shown that [54] 



- — =?=n 2 ' 

l—Ti 



With this result eq (90) can be combined with eq 
(89) to give 



i?=exp [-2n(kl/S) 1/ 



(91) 



This gives a simple relationship between reflectance 
and the fundamental optical parameters which 
has been shown to be valid for weakly absorbing 
samples [54]. 

It was found that the above simple particle theory 
disagrees with the Melamed theory for large refrac- 
tive indices [55]. The discrepancy has been attributed 
to a breakdown of the assumption that the externally 
incident radiation reflected at the particle surface 
is scattered equally in the upward and downward 
directions, an assumption which obviously becomes 
more questionable as n increases. By correcting 
the simplified particle model for this non-isotropic 
reflection of incident radiation, the relationship 



R _ l+r e (l-ft)/2-V(2a-a 2 )(l-7 e 74) 
r e /2+l — a 

was derived, which gives reflectances which agree 
well with those of the Melamed theory for large 
refractive indices [55]. 

Johnson [56] has pointed out that the Melamed 
theory [46] predicts that the diffuse reflectance 
decreases with increasing values of the relative 
refractive index, n. The same applies to the simplified 
particle model described above. Such behavior 
is reasonable for large values of n, but for values 
near unity it has been shown [56] experimentally 
that R— >0 as n— »1, as would be expected. 

The failure of the Melamed theory and the 
simplified particle model theory of Simmons for 
n values near unity is thought to be due to the 
invalid assumption that radiation is returned 
equally in all directions from a given particle [57]. 
It is found that for a rough-surfaced particle; i.e., 
one whose surface obeys the Lambert cosine law, 
three-fourths of the incident radiation reverses 
its direction while one-fourth does not [55]. It is 
assumed that the radiation which is incident at 
a particle surface may be divided into that which 
is randomly scattered and that which is transmitted. 
If that part which is scattered becomes isotropically 
distributed and the direction of the transmitted 
radiation is essentially unchanged, then the following 
expressions can be derived : 



r=^r e + {\- 



1 - 



-ft) (1 



-»G0 






=Jf.+(l-70(l-r.)Q«) 



1 



(i-HO-H' 



1— r t t 



1+2 7,1 



(93) 



(94) 



These equations for r and t are used with eq (43) to 
calculate the reflectance. It is found that the agree- 
ment with Johnson's data [56] is only qualitative, 
although the solution does have the feature that 
R — > as n — > 1. 



IV. Discussion of the Various Theories 

If we assume isotropic radiation and the simplest 
Gaussian quadrature n=l in eq (32), we obtain the 
Kubelka-Munk equations (16) and (17), with 



S= 



(j-\-a 



a , Tjr 2ot 

and K=- 



(T-\-OL 



Thus the Kubelka-Munk function may be regarded 
as a first approximation to the complete solution for 
the equation of radiative transfer. Detailed com- 
parisons show that the two solutions differ by no more 



than a few percent, and even this difference is of 
little consequence in a comparative- type measure- 
ment [12]. 

The Kubelka-Munk theory is in fact quite general, 
and encompasses many other two-constant theories 
which have been derived to suit certain select ex- 
perimental conditions [7,23]. With the Kubelka- 
Munk solution expressed in terms of hyperbolic func- 
tions, it is possible to write formulas for many specific 
applications, and to show the connection with several 
other theories. It was shown explicitly by Kubelka 
[21] that the various equations of Gurevic [20] and 
Judd [58] could be derived from the Kubelka-Munk 
equations. The Gurevic and Judd theories are less 
general, however, and do not encompass completely 
the Kubelka-Munk results. 

The Gurevic layer model is sufficiently general to 
allow us to use it in making some general observa- 
tions concerning the continuum models and their 
range of applicability. Since Kubelka [21] has de- 
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rived eq (12) of Gurevic in terms of the hyperbolic 
solution, we have a connection between the para- 
meters used in the two cases. In particular, 

L=(a 2 -l) 1/2 S=bS 

where a is defined by eq (18), and R m has the same 
meaning in the two theories. L and R m are related 
in turn to the Gurevic constants K x and K 2 through 
eqs (14) and (15). 

It will be recalled that K 1 and K 2 are positive 
definite constants which are characteristic of the 
light scattering medium (see eqs (10) and (11)). 
That is, we assumed thin homogeneous layers where 
T and R depend linearly on x. 

We find, in fact, that eqs (10) and (11), and cor- 
respondingly all continuum model results, are appli- 
cable only within a limited range. From eq (11) we 
have the obvious restriction that K 2 dx<^\. It is 
equally obvious that the equations do not apply 
when dx is made arbitrarily small, since the assump- 
tion of homogeneous layers requires it to be large 
compared with the dimension I of the scattering 
particles [14, 39]. Thus we have [16] 



l<<£dx<C 



K 2 



which can alternatively be written 
K 2 l<^\. 

We see that the attenuation due to absorption 
within a layer of thickness I must be very small 
compared with unity, and we are restricted to 
weakly absorbing and weakly scattering materials. 
It can in fact be shown by expanding eqs (12) and 
(13) in series as a function of x that eqs (10) and 
(11) follow only if one assumes that the medium 
scatters and absorbs weakly [59]. 

Several authors have investigated the relationship 
between continuum and statistical theories, ter 
Vrugt [60] compared the Kubelka-Munk and Bodo 
theories and found that the absorption constants 
determined by the two methods agree quite well. It 
was shown that for weak absorbers the parameters 
are related by 

Kx=kl (95) 

and 

Sx=^— (96) 

1— r v J 

If the crystals are not too irregular, the mean path 
length of radiation in the particle will be approxi- 
mately equal to the particle diameter and the 
absorption constants described by the two methods 
will be nearly equal. A more general relationship is 



exp (2kl)= 



2(1— r) 



F(R m )/F(r) 



(97) 



Poole [61] has shown that the Kubelka-Munk and 
Melamed theories lead to the same results for small 
values of the refractive index, n. In general, the 
absorption constants determined by the two methods 
are not the same and a curve is given by Poole 
which relates them as a function of n. It was sub- 
sequently shown by Karvaly [47] that the propor- 
tionality factor can be expressed analytically as 



2 A—r e 
3 n - 



l + r e 



(98) 



which is valid for materials of medium absorption 
as well. 



and that this relationship can in fact be used to 
estimate n for a powdered material. 

Several of the theories were compared by Com- 
panion [62] in an attempt to interpret reflectance 
spectra of metal oxides such as NiO and V 2 5 . It was 
found that the Kubelka-Munk theory gives unex- 
plained distortions and peak shifts, the Melamed 
theory requires a shading factor ^ ^0.1 rather than 
0.284 which is appropriate for spherical particles, and 
the Johnson estimate of the factor y is not a good 
approximation. These discrepancies are all thought 
to be a result of the failure of the theories to ex- 
plicitly include anomalous dispersion effects. 

Karvaly [47] has shown that the Johnson model is 
a limiting case of the Bodo theory which should be 
valid in the limit of large particles or materials of 
large optical thickness. Antonov-Romanovsky's for- 
mulae for regularly shaped particles are likewise 
derivable from the Bodo theory [47]. Careful meas- 
urements of the reflectance of two didymium glass 
powder samples (1=7.3 and 12.2 /zm) by Karvaly and 
Pinter [63] have shown that the Bodo and Antonov- 
Romanovsky theories give absorption constants 
which are accurate to about ±35 percent, whereas 
the Johnson and Melamed theories did not lead to 
satisfactory results in this case. 

As mentioned above, the Melamed theory fails to 
predict that R— >0 as n— *1, as must be the case since 
there are then no scattering centers in the medium 
[56]. In the Johnson and Bodo models the path of 
the radiation is altered by internal and external 
reflections and refractions at particle surfaces, where- 
as in the Melamed model it is diverted by "scat- 
tering" according to the Lambert cosine law. In all 
three cases it is assumed that Z^>\. This allows 
scattering to be neglected in the Johnson and Bodo 
models. The dependence of R on n in the Melamed 
theory is apparently due to the scattering of radiation 
internally incident on the particle surface which is 
taken to be independent of n. It is interesting to 
note that the three theories give about the same 
reflectance for refractive indices 1.7<^n<^2.0. 

It is of course not necessary to assume that the 
particles are rough-surfaced, i.e., that they scatter 
according to the Lambert cosine law. In the other 
limit they may be approximated as smooth-surfaced 
spheres whose reflection and refraction is governed 
by Fresnel's laws. This case has been treated by 
Simmons [64] and represents the opposite extreme 
between which most real systems lie. 
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The recently published modified particle model 
theory of Simmons [57] (Equations (93) and (94)) 
removes the discrepancy in the Aielamed treatment. 
This modified particle model theory is thought by 
Simmons to be the most nearly correct of all diffuse 
reflectance theories [39]. It ma}^ very well be true, 
but its general acceptance must of course await a 
thorough testing under a variety of conditions by 
various workers. 

The Fassler and Stodolski theory also shows con- 
siderable promise, but has not yet been thoroughly 
tested. It does lend itself to the investigation of 
important effects such as the influence of particle 
size distribution on reflectance properties [52, 65], 
and will probably receive considerable attention in 
the future. 



V. Interpretation of Reflectance Parameters 

One of the difficulties with weighing the relative 
merits of the various reflectance theories is the neces- 
sity of comparing them over a wide range of experi- 
mental conditions. In general, the data are rather 
incomplete and we must content ourselves with 
pointing out the physical significance of the param- 
eters used. We will do this in the present section, 
together with making some remarks concerning the 
characteristics of appropriate model systems. 

The internal transmittance of a powder layer in 
terms of hyperbolic functions has been given by 
Kubelka, [21] 



T<= 



(a sinh bSX+b cosh bSX) 



where a and b are previously defined (eqs (18) and 
(21)). If we assume that S is small compared with K, 



and 



A K 

T t = exp (-KX). 



Under these conditions the transmittance follows a 
Lambert-t}^pe law and K can be regarded as an 
absorptivity characteristic of the substance. 

The Kubelka-Munk constant K for a scattering 
medium is not equivalent to that which would be 
determined by transmission measurements on an 
identical material without scatter, however. For a 
beam traveling through the infinitesimal layer dx at 
an angle 6, the path length in the layer is dx/cosd, 
and thus the mean path length of radiation traveling 
downward is 

Jo I + ddcosd Uyj 



dx 



where dl+/dd is the angular distribution of the in- 
tensity in the positive x direction. We of course have 
an analogous expression for the /_ component. If 
we assume that the medium is an ideal diffuser the 
intensity is the same in all directions, and the angular 



distribution of intensity through a given plane 
becomes 

|±=J + sin 29 

so that eq (99) becomes equal to 2dx. Thus the ef- 
fective path length in the scattering layer is twice the 
normal la}^er width because of the random angular 
distribution, and the apparent absorptivity should be 
larger by a factor of two. 

The above result has been tested using colored 
filter glass by Kortiim and Oelkrug [66]. By measur- 
ing both the reflectance and scattered transmission 
of thin layers, both the scattering and absorption 
coefficients can be determined. Even though the 
powder cannot be packed to the same density as the 
original glass, the ratio of reflectance to transmit- 
tance values exceeded the theoretical limit, being 
in the range K/k~2.6-2.9. This difference was 
attributed to effects of total internal reflection which 
should increase the path length even further [8]. 
There are numerous factors which can cause devia- 
tions from the factor of two in non-ideal samples, and 
these factors have been discussed by Van den 
Akker [67]. We will continue to restrict our discus- 
sion to ideal systems, however. 

The statistical theories predict this ratio between 
absorptivity determined by reflectance and trans- 
mittance with varying degrees of success. As already 
mentioned, ter Vrugt [60] found them to be equal 
for weak absorbers (cf. eq (95)). The Antonov- 
Romanovsky theory [22] gives 

K = 1-Tq 
k l+r 

For the BG 24 filter glass used by Kortiim and 
Oelkrug [66] the ratio is 0.92, which does not agree 
well with the experimental values. 

For low refractive indices Poole [61] showed that 
the Kubelka-Munk and Melamed theories agree with 
a proportionality factor of 2.7, but much higher 
values are required for higher n values (cf. eq (98)). 
Fassler and Stodolski [68] have also derived a rela- 
tionship between K and k which involves the re- 
fractive index of the medium: 



1 



rtf 



k Su (1 — r e ) 

The ratio is calculated to be in the range 2.9-3.1 for 
the filter glass used by Kortiim and Oelkrug [66], 
which is in excellent agreement with the experimental 
values. The ^-dependence of K/k has also been cal- 
culated by Simmons [71] from his particle model 
theory. The result is 

K = 2n 2 
k 3 

which also agrees fairly well with the Kortiim 
and Oelkrug data [68]. 
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Multiple scattering in dense media is a phenomenon 
which is not completely understood. Mie [3] has 
given an exhaustive treatment of single scatter and 
has shown that the scattering varies from the A~ 4 - 
dependence of Rayleigh scatter for very small par- 
ticles to a X-independence for particles which are 
large compared with the wavelength of the radiation. 
There is no a priori reason to expect that these re- 
sults would apply quantitatively in the case of mul- 
tiple scatter, but it has in fact been found that 
there is a rather close correlation. Kortiim and 
Oelkrug [66] found that the scattering varies as 
^A~ 3 - 5 for 1<C\ as ^\ _1 for Z~\, and as ~\° for 
l^>\, where I is the mean particle diameter. It is safe 
to assume that the scattering coefficient S is inde- 
pendent of wavelength in most studies, but devia- 
tions are sometimes observed in the short wavelength 



end of the spectrum which are probably due to a 
breakdown of this assumption [69]. 

The scattering coefficient has been found in several 
studies to be inversely proportional to the mean 
particle diameter I, [70-72] the proportionality factor 
depending apparently on the nature of the material 
being studied. Using the simplified particle model, 
Simmons was able to derive the result [69]. 

This should not be taken as much more than an order 
of magnitude estimate, since a proportionality factor 
somewhat larger than unity is apparently more ap- 
propriate in many cases [66]. The modified particle 
model gives [39] 



S=l~ 



{ 



7 °+\ (l-rO (1-70 \j±r 



(l-r«/2) (1-7./2) 



l+ftt/2 



]}■ 



Although this result appears to not yet be thoroughly 
tested, the term in braces is said to be relatively 
insensitive to changes in absorptivity [41], and 
probably varies significantly only with changes in 
refractive index, n. The same is true of the Antonov- 
Romanovsky result [22], 



S-- 



vi W 



and that derived by ter Vrugt [60] (see eq (96)). 
Quantitative calculations with the Antonov-Roma- 
novsky expression give scattering constants that are 
too small by about a factor of three [47]. 

By taking the particle size distribution into ac- 
count, Fassler and Stodolski [65] were able to cal- 
culate scattering coefficients which agree very well 
with those determined experimentally. Their ex- 
pression is 



Sl'-- 



l+r*-g 
l-r*-£~ 



-a-(2f-l)(l-a) 



l+r*+g 



-(2/-l)(l-a) 



(100) 



where a and f t have their previously defined mean- 
ings, it being recognized that the factor (2/-1) 
defines the anisotropy of radiation in the particle 
layer. £ is the effective cross section for holes in the 
particle layer, r* is the mean external reflection 
coefficient for side scatter taking account of holes 
and back scatter, and ~V is a parameter which is 
related to the mean particle diameter. For a particle 
size distribution function of the form 



V is given by 



f(l)=Al B e- Bl/1 



7/- g±8 7 

1 ~ B 



Fassler and Stodolski [65] emphasize the fact that 
the scattering is particle size-dependent, so that for 
quantitative reflectance spectroscopy it is essential 
to have a set of standards with known, small particle 
size distribution and of various mean particle di- 
ameters. Then it is possible to mix the sample with a 
standard of similar particle size so that the scattering 
characteristics are not greatly changed. This is 
important for comparative measurements. 

The effect of particle shape is also important. It 
appears at the present time that there is no better 
way to test the various reflectance theories than 
to make transmission and reflectance measurements 
on a given glass sample before and after grinding, 
respectively. The characteristics of glass samples 
subjected to various grinding conditions have been 
studied by Karvaly and Pinter [73]. These factors 
will not be discussed here, but it should be em- 
phasized that the proper definition of reflectance 
standards must certainly quantify not only size and 
shape, but the precise nature of any adsorbent 
materials and the surface conditions as well. 
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